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Abstract 

Following the general scheme presented in [Am. J. Phys. 78 (8), 842-850 (2010)] we analyze the 
quantum-mechanical behavior of a system described by a one-dimensional asymmetric potential 
constituted by a step plus i) a linear barrier or ii) an exponential barrier. We explicitly solve the 
energy eigenvalue equation by means of the integral representation method, which allows us to 
classify the independent solutions as equivalence classes of homotopic paths in the complex plane. 

We discuss the structure of the bound states as function of the height Uq of the step and we 
study the propagation of a sharp-peaked wave packet reflected by the barrier. For both the linear 
and the exponential barrier we provide an explicit formula for the delay time t{E) as a function of 
the peak energy E. We display the resonant behavior of t{E) at energies close to Uq. By analyzing 
the asymptotic behavior for large energies of the eigenfunctions of the continuous spectrum we also 
show that, as expected, t{E) approaches the classical value for E — oo, thus diverging for the 
step-linear case and vanishing for the step-exponential one. 

PACS numbers: 02.30.Uu, 02.30. Gp, 02.30.Hq, 03.65. Ge, 03.65.Nk 
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I. INTRODUCTION 



In a previous paper^ we have analyzed the quantum-mechanical behavior of a system 
described by a one-dimensional asymmetric potential formed by a step plus a harmonic 
barrier (the "step-harmonic" potential), by using the integral representation method.- We 
investigated the behavior of the discrete energy levels (as a function of the height of the 
step) and of the delay time r of a wave packet coming from infinity and bouncing back on 
the harmonic barrier, as a function of the packet's peak energy and of the height Uq of the 
step. 

Among the convex or concave locally bounded symmetric and confining potentials the 
harmonic oscillator is a threshold, in that it gives rise to classical isochronous oscillations 
and thus to evenly spaced energy levels when quantized.-*^ 

In our quantum mechanical step variant of the problem we recover both these features 
in the limit in which the step height Uo oo, and the potential reduces to the half-space 
harmonic oscillator. Then, it is conceivable that the harmonic one is the only confining 
barrier which displays a constant nonvanishing delay r in the limit of high energies. For 
steeper barriers we expect r to vanish at high energies, while for milder ones we expect 
the delay to become infinite in this limit, in accordance with the corresponding classical 
situations. Similarly, we expect that, as f/o — C)0, the spacing between two neighboring 
discrete levels tends to infinity in the former case and to zero in the latter. 

In the present paper we corroborate this conjecture by analyzing two examples, the "step- 
linear" potential and the "step-exponential" one. Both these problems can be exactly solved 
by the use of the integral representation method. 

II. THE STEP-LINEAR POTENTIAL 

Consider the "step-linear" potential 




Mx X <0 



(1) 
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Fig. 1: The "step-linear" potential. 

where M and Uq are positive parameters. If E denotes the energy of the particle and m its 
mass, the time-independent Schrodinger equation is 

cPuix) 



2m dx^ 



+ U{x)u{x) = Eu{x) 



This, for a; < 0, writes 



d'^u{x) /2mM 2mE\ 

dx'^ \ h? J 



u{x) = . 



It is convenient to define 



a :- 



/2mM\ 



1/3 



which allows us to recast Eq. (E]) as follows: 

d^u{y) 



(3 



aE 



y := ax 



+ {y + /3)u{y) = 



which is the Airy equation.-i^ The general solution of Eq. ([5]) is 

u{y) = C Ai{-y - f3) + D Bi{-y - (3) 



(2) 



(3) 



(4) 



(5) 



(6) 



where C and D are arbitrary integration constants, and the two linearly independent solu- 
tions Ai and Bi are expressed in appendix |X] in terms of the integral representation method. 

Since Bi(x) — )■ +oo for x — )■ +oo, in order for Eq. to be an eigenfunction, we must set 
D = 0. Hence, we obtain 

u{x) = CAi{-ax - P) . (7) 



For X > 0, the solution of Eq. ([2]) has the following from: 

{^^ikx Q^-ikx E>Uo 
(8) 
Fe"^^ E <Uo 

where A, B and F are arbitrary integration constants and 

hk := ^2m\E-UQ\ . (9) 

A. The case E < Uq: bound states and level spacing 

If E < Uo we obtain 

{CAi(-ax - (5) a; < 
(10) 
P^-kx ^ ^ Q 

The requirement of continuity of u{x) and of its first derivative in x = is expressed by 

CAi(-/3) - F = 

(11) 

Caki{-(5) -Fk = 
System (ITTl) has a non trivial solution iff 

^^V^, (12, 

where /3o := aU^/M. The energy levels are determined graphically by the intersections of 
the curves at the two sides of Eq. f[T2]) . An example is depicted in Fig. [2] for /3o = 6. In 
the limit f/o — t- oo the step of Eq. ([1]) becomes an infinite barrier. In this case, the energy 
levels correspond to the zeros /3„ of Ai(— /3), the denominator of Eq. ( IT2|) . As expected, these 
energy levels are the ones of the symmetric confining potential U{x) = M\x\ corresponding 
to the odd eigenf unctions of the latter (appendix iBl). 

It is instructive to study how their separation behaves for large energies. To this purpose, 
consider the asymptotic expansion of Ai(— /3) for large values of 



Ai(-/3) 



sm 



oo oo 

+ 9 - COS (c + 9 $:(-)'=c,.+iC"^'^-^ 



k=0 k=0 



(13) 
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Fig. 2: (color online) Graphical solutions of Eq. (jl2p for the energy levels. The red solid line 
represents the graph of the function Ai'(— /3)/Ai(— /3) while the blue dash-dotted line represents 



the graph of -y/^^o^^. In this plot, we have chosen (3q = 6. 

where ( := 2/3^/^/3 and the coefficients of the series expansions are given by 

r(3t + i) 

' 54»A!r(*:+i) ■ ^ ' 

Note that expansion ( fT3i) can be proved by easy but rather tedious calculations starting 
from Eq. (|A5l) . 

The inversion of the asymptotic expansion (fT3l) allows us to find, for large values of (3, 
the following approximate solution of the equation Ai(— /3) = :- 

= T(t„) , with = Iti (4n - 1) , (15) 



T{x) ~ x'/' ( 1 + __-+...) , X -> oo . (16) 



and 

5 1 5 1 
48x2 ~ 36x 

Thus, at the leading order of Eq. f|T6l) the approximate zeros have the form 

2/3 



-vr {An - 1) 



(17) 



This is an excellent approximation to the zeros of Ai(— x), already from the first one (see 
Table H]). For n — )■ oo we get from Eq. ( IT7|) . to leading order in 1/n, 



n 


/3„ (exact) 


/3„ (approximate) 


Relative Error 


1 


2.33811 


2.32025 


0.76 X 10-2 


2 


4.08794 


4.08181 


0.15 X 10-2 


3 


5.52055 


5.51716 


0.62 X 10-3 



TABLE I: The first three zeros of the Airy function (n = 1, 2, 3) compared with the corresponding 
approximate zeros from Eq. (jl7|) . 



1/3 

3n 



n oo 



(18) 



Therefore, the level spacing decreases with n and vanishes in the limit. The spacing behavior 
n~^^^ is the threshold between concave and convex potentials. 



B. The case E > Uq: scattering and delay 

In this case, the (improper) eigenfunctions have the form 



uix] 



CAi{-ax - f3) x<0 



(19) 



(20) 



The junction conditions in x = are: 

CAi(-/3) =A + B 
-CaAi{-f3) = tk{A- B) 

Solving for the constants, the normalized (with respect to k) improper eigenfunctions are 
given by: 

1 f n[/3(fc)]Ai[-ax - (3{k)] x<0 

y 2,71 I ^—ikx _|_ ^ikx+iS(k) X > 

where 

-1 iA;Ai(-/3) - aAi'(-/3) 



(21) 



a 



n(/3) = 2 Ai(-/3) + T^Ai'(-/3) 



ikki{-(5) + aAi'(-/3) 
As expected, the continuous part of the spectrum {E > Uq) is simple. Note that 

'aAi'[-f3{k)]' 



5{k) = 2 arctan 



k Ai[-/3(A;)] _ 



(22) 



(23) 



-j:E{k)t 



From Eq. (12T]) a generic wave packet 

ip{x,t) = I dk c{k)uk{x)e^^^ 
Jo 

has the form 

1 \ dk c{k)U[f3{k)]Ai[-ax - f3{k)]e~T^^^''^' x<0 
ip{x, t) = —= < 

[ /o" dk c{k) [e'kx+:5{k) ^ g-ifcxj = ^^^^ + ^.^ a; > 

Then, writing c(A;) = |c(fc)|e*'^*-'^'', ■j/'m and ■j/'refi take the following form: 



(24) 



(25) 



^i„(x,t) 
V'refl(a;,t) 



+ 00 



/27r JO 



dk |c(A;)|e-*['=^+^('=)*-^('=)l , 



where 



(26) 
(27) 

. . (28) 
If c{k) is sufficiently regular and non-vanishing only in a small neighborhood of some k, 
then ipi^ and ^/'refl represent wave packets which move according to the following equations 
of motion:-!^ 



dVL 
~dk 



t + 



k=k 



d'-f 

dk 



hk . . p . , 



k=k 



m 



m 



for the "incoming" wave packet, and 

Xve& 



dQ 


t - 

k=k 


d^ 




d6 


- P. 


dk 


dk 


k=k 


dk 


k=k m 



, . m d5 
p dk 



k=kJ 



(29) 



(30) 



for the reflected "outgoing" one. 

The solution thus built represents a particle of well defined momentum p = hk which 
approaches the origin from the right, interacts with the linear potential (at t = to)^^? and 
is totally refiected. The phase shift results in a delay r in the rebound,J^ caused by the 
interaction with the confining linear barrier. From Eqs. (jl]) and (jH]) it follows that 



. ~, ah d6 



(31) 



where (3 := (3{k). We compute r from Eq. (1231) . Taking into account the Airy equation 
Ai"(-/3) = -/3Ai(-/3), we obtain 



ah 



1 , . Ai(-/3) Ai-(-/j) 
2(/3-/3o)^^Ai'(-/3)^ Ai(-/3) 



'Ai'(-/3) 



Ai(-/3) 



(32) 



where we have suppressed the tilde on the packet peak energy (3. 

We are interested in the behavior of the interaction time for large values of /3, i.e. for 
incoming packets with high energy. To this end, we need the asymptotic expansion of 
Ai'(-/3) for 13 ^+00 



Ai'(-/3) = - " 



cos 



c 



k=0 ' k=0 

where ( := 2/3'^''^/3 and the coefficients dk are 



(33) 



with Ck given in Eq. fll4p . 

Then, dividing the two asymptotic expansions of Ai'(— /3) and Ai(— /3) we obtain to leading 
order in /3 

Thus, Eq. (132|) becomes 

r(/3)~2— v^, /3->oo. (36) 

Hence, reintroducing the physical variables, the high-energy behavior of the interaction time 
is 

r{E) ~ , (37) 

which is exactly the time a classical particle arriving from infinity with energy E would spend 
in the a; < region. In Fig. [3] we plot t(/3) for a choice of different values of Po- Note the 
resonances located at the points /3 ~ ?7„ (n = 1, 2, . . . ), zeros of Ai'(— /3), corresponding to the 
formation of metastable states at the respective energies En — Mrjn/a. The values Mrjn/a 
are the energies of the excited states of the confining potential M\x\ (see appendix [B|) . 
corresponding to even eigenfunctions. The resonances have lifetimes which decrease as 
the corresponding energies increase and move farther away from the threshold energy Uq. 
Conversely, as Uq increases, the lifetime of the resonance closest to the height of the step 
becomes progressively longer and then infinite when the resonance turns into the next bound 
state. This behavior is evident in Fig. 121 in which the first three plots correspond to values 
of Pq for which there is only one bound state. In the successive three plots the resonance at 
/3 ~ r/i has disappeared, having turned into the second bound state. 
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Fig. 3: (color online) Plots of the delay r (in units of ha/M) versus the energy /3 of the incoming 
wave packet (red solid lines) for six values of the step height: /3o = 1-5, 2.1, 2.7 (upper panels) and 
/3o = 3.5,4.0,4.5 (lower panels). The parabolic blue dotted lines represent the classical delay, given 
in Eq. (l37l). 



Comparing Fig. |2]with figure 6 of paper 1| we note that, whereas in the step-harmonic 
case the graph of t(/3) oscillates with decreasing amplitude about the straight line t = n/u 
(the half period of the oscillator), in the step-linear case the corresponding graph similarly 



oscillates about the parabolic line t{E) = 2y/2mE/M, corresponding to the delay of the 
classical particle. Furthermore, whereas in the step-harmonic case the resonances are evenly 
spaced, in the step-linear case their spacing decreases with the energy, corresponding to the 
behavior as a function of the energy of the eigenvalues of the corresponding (symmetric) 
potentials U{x) = and U{x) = M\x\. 



III. THE STEP-EXPONENTIAL POTENTIAL 

Consider the "step-exponential" potential 



U{x) 



K {e-'-"/" - 1) X < 



(38) 



X > 
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Fig. 4: The "step-exponential" potential. 

where k, a and Uq are positive parameters. For x < 0, it is convenient to introduce the 
following dimensionless quantities: 



irana 



a 



in terms of which the time-independent Schrodinger equation writes as 

2(Pu{z) du{z) 



z 



+ z^ + (/3 - z') u{z) = 



Setting := —f3, Eq. f HOl) can be cast in the form of a modified Bessel equation:-"^ 

^d'^ui^z) ^ du{z) 



dz^ 



+ z- 



dz 



+ z^) u{z) = 



whose general solution is^^ 



u{z) = C K,{z) + D h{z) , 



(39) 



(40) 



(41) 



(42) 



where C and D are arbitrary integration constants and and Ip are the modified Bessel 
functions of order u = i\f^. 

The function lui^z) diverges exponentially for z — ?■ +oo.^ For this reason, in order for uix) 
to be a proper (or improper) eigenfunction, we must set D = 0. Therefore Eq. ( 1i2|l reduces 
to 

u{x) = C K,^ [ae-^l^") . (43) 

Like the Airy functions, the solutions of Eq. fHT]) can be conveniently studied by means of 
the integral representation method (see appendix [C] for details). 
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Fig. 5: (color online) Graphical solutions of Eq. ()46|) . The red solid line is the graph of 
K'.^{a) / K-^{a), whereas the blue dash-dotted line represents \/ I3q — f3/a for the choice /3o = 24 
and a = \. 



A. The case E < Uq: bound states 

If E < Uq we obtain 

u{x) = I ' ~ . (44) 

The junction conditions in x = give 

= '^^kK.^ia) , (45) 
which, setting (3q := 8m{Uo + K)a^/h'^ can be recast in the form: 

^^v^(") ^ (46) 

Graphical solutions of Eq. (H6!) are shown in Fig. |5l Analogously to what happens in the 
step-linear case [Eq. f lT2|) ]. in the limit of an infinite barrier {Uq — ?■ oo) the energy levels are 
specified by the zeros of K^^{a), the denominator of Eq. (H5|l . as a function of /3 and they 

11 



are the ones of the symmetric confining potential U{x) = k (el^l/"" — l), corresponding to 
the odd eigenfunctions of the latter (appendix [Dll . 

To study how the energy levels behave for large energies, we employ the following formula 
for the asymptotic behavior of the function K^^{a) for large 

„2 



2tx 
7^ 



a 



4v^ 



- V/3+ V/31o, 



TT 

a ) A 



1 + 



(47) 



for /3 — i- oo . Note that expansion (H7|) can be proved starting from Eq. (ICSp . Therefore, 
the zeros of K^^^a), as a function of /3, are asymptotically the solutions of the following 
equation: 



a 



rnr . 



Solving for /3„ we obtain 



(3n 



exp 



2W 



2mr 



ae 



(48) 



(49) 



where W{x) is the Lambert function.-. Then, since W{x) ~ logx — log log x for x — )■ oo, 
from Eq. fHHj) we have for large n that 



/3. 



(50) 



We see from Eq. (|50l) that the potential U{x) = k (e'^'/°" — l) behaves for large x as an 
infinite square well whose width, up to inessential factors, grows as logn, an intuitive fact. 
Moreover, 

for n — )■ cxD, proving thus that the level spacing diverges. 



B. The case E > Uq: scattering and delay 



The unbound eigenstates have the form 



CK. miae-'^l'^") a; < 
u{x)={ ' - , (52) 



and the junction conditions are 



(53) 

CaK'.^{a) = i2ak {B - A) 
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Therefore, the normahzed (with respect to k) improper eigenfunctions are given by: 



Uk(X 



27r 



+ e 



ae 



where hk = ^J2m{E — Uq) and 



Hence, 



Oi 



a] 



-x/2a^ X < 

X > 



2ikaKi^{a) 



(54) 



2ikaKi^{a) + aK'. 



a] 



6{k) = 2 arctan 



K' , (a) 



(55) 



(56) 



Then, following the same argument adopted for the step-linear case, we obtain the following 
formula for the delay r of the rebound of an incoming wavepacket with peak energy /3: 



h d(3 



1 



d K' (a) 
+ -j^ log ■ 



2(/3-/3o) "i^.v?(«) 



Using Eq. ( H7|) . we obtain for large values of f3 



a] 



-^cot(-v^+V^log^ + ^ 



from which 



d/3 Ki^(a) 



1 

2^ 



1 



sm 



1 , 2v^ 
log 



a 



(57) 



(58) 



(59) 



2(-v^+y^log^ + f^ 
A comment is here in order. In general, taking the derivative of an asymptotic expansion 
with respect to the variable or a parameter may lead to wrong results. However, in our case 
this procedure can be justified using the integral representation of Eq. (\C9\} (we leave this 
as an exercise for the interested reader). 

Thus, inserting Eq. ( 13^ into Eq. we obtain the asymptotic behavior of the delay 
time for large /3's, namely 



6ma 



log 



2v^ 



a 



or, in terms of the energy of the particle 

t{E) - 



2(T \og{2E/K) 



(60) 



(61) 



^j2E/m 

As expected, Eq. fl6T]) coincides with the large energy value of the half period of the classical 
particle subjected to the confining potential U{x) = k {e^^^/'^ — l) . 
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Fig. 6: (color online) Plots of the delay r (in units of 8ma'^ /h) versus the energy (3 of the incoming 
wave packet (solid red lines) for six values of the step height: (3q = 10, 11, 12 (upper panels) 
/3o = 22,24,25 (lower panels). The blue dotted lines represent the classical delay. 

IV. CONCLUSIONS 



Regarding the structure of the discrete energy spectrum as a function of the height Uq 
of the barrier, in the two examples treated in this paper [Eqs. (Tj) and (|38l) ]. the same 
considerations apply as those of the concluding section of paper The only difference 
is that the energy levels fkj{n + 1/2), n = 0, 1,2,3, of the harmonic oscillator have to 
be replaced here by the corresponding levels En of the confining potentials f lBl|) and (IDip . 
respectively. In the case of potential ([1]) the level spacing goes to zero as the energy increases 
[Eq. (|T8l) ]. while in the case of potential (|38l) it approaches infinity [Eq. (ISTl) ]. As regards 
the continuous spectrum, we provide in both examples exact expressions for the delay of a 
wavepacket reflected from the barrier, as a function of the peak packet energy [Eqs. (!32l) and 
( l57|l ]. As expected, in both examples these delays exhibit a series of resonances for energies 
not much larger than Uq, while for large energies, they approach the classical values. The 
step-harmonic potential is a threshold separating the potential barriers for which the delay 
time goes to inflnity at large energies from those for which it vanishes. 
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An entirely similar discussion can be applied to the step variant 




(62) 



of any symmetric potential V{x) (y{x) = V{—x)) such that lim^^-too ^(2;) = +00. Indeed 
the energy eigenvalue equation for V{x) has two linearly independent solutions ul{x) and 
vl{x) the first of which approaches zero very rapidly as x — > —00 whereas the second one 
diverges steadily without oscillating, and two linearly independent solutions ur{x) and vr{x) 
having a corresponding behavior for x — )■ +00.-*^^ Since ul{x) = a{E)uji{x) + b{E)vR{x), 
the energy eigenvalues are the roots En of the equation b{E) = 0. Since the potential is 
symmetric, these roots correspond to even and odd eigenfunctions alternatively, the ground 
state being even. However, in the general case the eigenvectors cannot be found explicitly. 
Therefore, for example, no explicit formula is available in general for the delay time of the 
reflected packet in the corresponding step variant potential ( l62l) . 
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Appendix A: Derivation of the Airy functions by the integral representation method 

We look for a solution of Eq. (j5]) of the form 



where 7 is a generic path in the complex plane and / is a function holomorphic in a suitable 
domain. Inserting Eq. (lAip into Eq. ([5]) we find 




7 



(Al) 




(A2) 



Integrating by parts, we obtain: 




7 



(A3) 



Therefore, E{y) is a solution of Eq. ([5]) if 








and 




(A4) 
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Fig. 7: (color online) Possible paths of integration, each one corresponding to a solution of Eq. 



Hence, a class of solutions of the Airy equation is of the form 



Eis{y) = dt exp 
J y 



3 



| + + 



(A5) 



where 7 is a suitable path for which the contour term vanishes. 

The integrand of Eq. flA5|) is an entire function; therefore, by virtue of Cauchy theorem, 
every closed path represents the trivial solution E{y) =0. 

Consider an unbounded path. In order for [e*^/(t)]g^ to vanish, we require the leading 
term in the exponent of f{t) (i.e. t^) to have a negative real part. Therefore, the acceptable 
unbounded paths are those whose phase is confined to the regions | < + ^rni < | 
{n = 0,1,2). These possible paths are showed in Fig. [71 where the allowed sectors | < 
(p + |n7r < I (n = 0, 1, 2) are shaded. 

Paths with both endpoints in the same sector (e.g. F4 in Fig. [7]) can be closed at infinity 
using Jordan's Lemma; therefore, they correspond to the trivial solution. The only non- 
trivial paths are those which link different sectors. There are only 3 non-equivalent classes of 
such paths which we dub Fi, F2 and F3 respectively (see Fig. [7j). Taking into account Cauchy 
theorem, these paths satisfy the relation Fi + F2 = F3 in the sense that the corresponding 
solutions are not independent. The conventional Airy functions Ai[z) and Bi{z) are the 
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independent solutions of w"{z) — zw{z) = such that^ 

3-2/3 



Ai(0) 
Bi(0) = VS 



r(2/3) 

3-2/3 



Ai'(O) = - 
Bi'(O) = VS 



3^ 

rws) ' 

3-1/3 



r(2/3) ' V ^ V Y{l/3) ■ 

Denoting by E'^^\y) the solutions in Eq. ( lASP corresponding to the paths Fj (i = 1, 2), it is 
not difficult to show that 



Bi(-y-/3) 



2tt 



Ef{y)-Ef{y) 



We leave the details^ to the interested reader {hint: Check the above expressions and their 
first derivatives in Eq. (lASp for y = [3 = 0. In this case the integrals Eq\o) correspond to 
Euler r functions). 



Appendix B: The confining symmetric linear potential 



Consider the confining symmetric potential (see Fig. |8]) 



U{x) = M\x\ . 



(Bl) 




Fig. 8: The "linear" well. 
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The eigenf unctions can be written as 

u(x) = CiAi(—ax — B) for a; < 

^ ^ ^ ^ , (B2) 

u{x) = C2Ai{ax ~ f3) for x > 

where Ci and C2 are integration constants to be determined by means of the junction 
conditions in x = 0. The latter can be resumed as follows: 

• If Ai(— /3) 7^ 0, then the requirement of continuity of u{x) in x = implies that 
Ci = C2. Moreover, the continuity of the derivative implies Ai'(— /3) = 0. This 
condition determines the even eigenfunctions and the corresponding eigenvalues. 

• If Ai(— /3) = 0, then the continuity of the derivative implies Ci = — C2. This condition 
determines the odd eigenfunctions and the corresponding eigenvalues. 

Appendix C: Derivation of the modified Bessel functions by the integral represen- 
tation method 

Consider the modified Bessel equation 

z\"{z) + zu\z) - (z/' + z^)u{z) = , (CI) 
with 2; > 0. A convenient kernel for the integral representation is the following: 

= ^''e-"™'^* . (C2) 

We look for solutions of the form 

u{z) = j dt f{t)K{z,t) . (C3) 



. (C4) 



Plugging Eq. (lC3j) into Eq. ((CTJ), we get 



z'^' / dt 



7 



/(t)^ (sinh(t)e-"™^'^*) + 2uf{t) cosh(t)e-"™^*^* 



Integrating by parts, Eq. (IC4I) gives 



[z"+^e-"^°^^V(t) sinh(t)]g^ - z"^^ / dt [f\t) sinh(t) - 2uf{t) cosh(t)] e"^™^^* = . (C5) 
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Im(t) 












27r 
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R.e(i) 




Fig. 9: (color online) Cut of the complex t-plane. The dashed red and the solid blue thick paths 
represent the only classes of paths which provide independent solutions to the modified Bessel 
equation. 



Then, the integral at the right hand side of Eq. (1C3I) is a solution of Eq. (ICll) if 



f{t) sinh(t) - 2z//(t) cosh(t) = 



and 



v+i ^g-.coshtj^^^ sinh(t)]^ = . (C6) 
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Up to a normalization, the solution is f{t) = sinh(t)^'^. If u is integer then f{t) is either 
entire {u > 0) or meromorphic {u < 0), otherwise it has infinite branch points located at 
tn = inn = 0, ±1, . . . ), see Fig. [91 In the latter instance, the usual procedure is to define 
a domain in which the function is holomorphic by cutting the t-plane and thus forbidding 
loops around the branch points. In Fig. [9] a convenient choice for the cuts is also shown. 
Recalling that Iie{z) > 0, the contour condition z'^'^^ [e~^™'^^''*^/(t) sinh(t)]g^ = writes 



(C7) 



There are 4 different classes of paths for which Eq. (1C7I) is satisfied and 

f smh{tfe-''°''^'z''dt , 

J y 



UiZ 



(a 



is well defined. The "paths zoology" is more complicated and rich than in the linear and in 
the harmonic case.- 
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1. Closed paths. For any closed path, the contour condition is trivially satisfied and any 
integral along a path enclosing a region where the integrand function is holomorphic 
(i.e. the path does not cross the cuts) vanishes. An example is shown by the thin 
dashed black line in Fig. O 

2. Infinite paths. For paths whose endpoints are both at infinity, the function 
sinh(t)^'^"'"^ diverges or oscillates. On the other hand, the exponential e"^^"*^*^* van- 
ishes for Re(cosht) — )■ +00 (recall that z > 0). Since Re[cosh(x + iy)] = cos(?/) cosh(x) 
then 7 must stretch at infinity in one of the sectors defined by —tt/2 + 2mr < lm(t) < 
7r/2 + 2r27r {n = 0, ±1, ...), which are represented by the shaded regions in Fig. [91 Inci- 
dentally, in these bands, when there are no cuts, one can "close" the paths at infinity, 
by virtue of Jordan's Lemma. Examples of this class of paths are the black solid thin 
lines of Fig. |9l 

3. Semi-infinite paths. By "semi-infinite" paths we mean paths starting from a point, 
say to, and ending at infinity. These paths must go to infinity in the shaded bands 
— 7r/2 + 2n'K < lm{t) < 7r/2 + 2r27r. For the starting point t^, the contour condition 
demands that sinh(to)^'^'''^ = 0. This means that these paths must start from one of 
the points t„ = mvr, which are the zeroes of the hyperbolic sine function. It is easy to 
prove that the integral of Eq. flC8|) performed along any two such paths lying in the 
same band gives the same result (indeed, recall that sinht is periodic). Two examples 
of this class of paths are the blue solid thick lines in Fig. |9l 

4. Finite paths. These are the paths starting and ending in two points say and tf, 
with < 00. The contour condition can be satisfied in two different ways: either 
the values of the contour part are equal at the endpoints, or the contour part vanishes 
at the endpoints. The former case accounts for paths which do not cross any cut and 
start from any point tj, ending at tf = ti + 2nm (n = 0, ±1, . . . ). Examples of this 
class of paths are represented by the red dash-dotted thick line in Fig. [91 The latter 
case is realized by paths connecting the branch points and is represented by the red 
dashed thick line in Fig. [91 

Taking into account the periodicity of the integrand function (the period is 27ii in the domain 
where it is holomorphic), and the Cauchy theorem, it is easy to show that the integrals along 
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the two kinds of finite paths are proportional to the integral along the "fundamental" path 
[Ojivr]. Moreover, the integrals along any one of the infinite paths are linear combinations 
of the ones performed along the finite and semi-infinite paths. 
In conclusion, two linear independent solutions of Eq. (ICip are 



where K,y and I,y correspond, respectively, to the integrals performed along the solid blue 
and the dashed red thick lines in Fig. |9] (a semi-infinite path and a finite one). It can be 
shown that both solutions (derived here for z > 0) can be analytically continued throughout 
the whole 2;-plane cut along the negative real axis.- 

Appendix D: The confining symmetric exponential potential 

Consider the confining symmetric potential (see Fig. [TO!) 




(C9) 



and 




(CIO) 



U{x) = n (el"!/'^ - 1) , 



(Dl) 



whose eigenf unctions can be written as 




U{x) 



X 



Fig. 10: The "exponential" well. 




(D2) 
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where Ci and C2 are integration constants to be determined by means of the junction 
conditions in a; = 0. The latter can be resumed as follows (compare the analogous situation 
for the linear well in appendix |B]) : 

• If K^^[a) 7^ 0, then the requirement of continuity of u{x) in x = implies that 
Ci = C2. Moreover, the continuity of the derivative implies K'.^{a) = 0. This 
condition determines the even eigenfunctions and the corresponding eigenvalues. 

• If K^^{a) = 0, then the continuity of the derivative implies Ci = — C2. This condition 
determines the odd eigenfunctions and the corresponding eigenvalues. 
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Note that the argument of the complex valued function c{k) determines to. 

The delay is calculated with respect to the case of instantaneous reflection, which takes place 

in presence of an infinite barrier and for which 6{k) = tt. 
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